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O ' Abstract 

■ A central problem in quantum computing is to identify computational tasks which can be 

^ ' solved substantially faster on a quantum computer than on any classical computer. By studying 

the hardest such tasks, known as BQP-complete problems, we deepen our understanding of the 
power and limitations of quantum computers. We present several BQP-complete problems, 
. including Phase Estimation Sampling and Local Hamiltonian Eigenvalue Sampling. Different 

than the previous known BQP complete problems (the Quadratically Signed Weight Enumer- 
ator problem and the Approximation of Jones Polynomials 12, fSlEDi our problems are of 
a basic linear algebra nature and are closely related to the well-known quantum algorithm and 



> 

. quantum complexity theories. 

^ . 

O ■ 1 Introduction 

, The last decade has witnessed a great surge of fruitful studies in the new paradigm of quantum 

computing, especially after Shor discovered the celebrated quantum algorithms for factoring 
Oh! and discrete logarithm |25| . While many other areas like quantum complexity theory, quan- 

tum cryptography and quantum error correction have been rapidly developed, the progress in 
^ I quantum algorithm design, one of the most central tasks of quantum computing, appears much 

slower than what people had expected. Especially, designing quantum algorithms that have 
super-polynomial speedup over their classical counterparts seems very hard, and we have only 
succeeded on few problems, including Hallgren's polynomial time algorithm for Pell's equation 
and class groups van Dam, Hallgren and Ip's polynomial time algorithm for some hidden 
^ ' shift problems [55], Kuperberg's subexponential time algorithm for HSP over the dihedral group 

H ; ESI, and some others. In a recent survey by Shor |27j . two major possible reasons are proposed 

to explain the difficulty. One is that we have no good intuitions to design quantum algorithms 
due to our lack of quantum experience. The other is rather pessimistic: there may actually be 
only a few problems for which quantum computers have significant advantages over classical 
computers, and probably we have already found almost all of them. 

Considering the enormous payoff of efficient quantum algorithms and these possible reasons 
for the difficulties to discover them, it is central to identify the class of problems that quantum 
computers can have substantial speedup over classical computers. Stated in the language of 
complexity theory, the task amounts to pinning down BQP, the class of languages decidable 
by a uniform family of polynomial-size quantum circuits with bounded error. While the class 
could be equivalently characterized by various models, including the quantum Turing machine 
0, quantum circuits |82j and quantum adiabatic computation [111 Bj. identifying it is difficult. 
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just as it is difRcult to characterize other important classes such as NP. Nevertheless, it 
was recently found that the Quadratically Signed Weight Enumerator problem 22 and the 
problem of approximately evaluating the Jones polynomial j^l I13L E] are BQP-complete^. 
Such problems are the hardest problems in BQP. Consequently, the study of BQP-complete 
problems is important as it helps us to deepen our understanding of BQP by fully capturing 
the power of quantum computation. This is similar to the study of NP-complete problem which 
aims at understanding the power of nondeterminism. 

In this paper, we provide several other BQP-complete problems, including the Local Hamil- 
tonian Eigenvalue Sampling (LHES) problem and the Phase Estimation Sampling (PES) prob- 
lem, which are (roughly) to sample an eigenvalue of a given Local Hamiltonian or a unitary 
matrix. Both problems are of an elementary linear algebra nature, which may make them acces- 
sible to more computer scientists. More importantly, these two problems are very closely related 
to well-known quantum algorithm and quantum complexity theories. First, the Local Hamilto- 
nian problem, which basically estimates the minimal eigenvalue of the given local Hamiltonian, 
is well known to be QMA-complete Our result on LHES says that if we do not 

aim at estimating the minimal eigenvalue, but at sampling an eigenvalue according to a very 
natural distribution, then the problem becomes BQP-complete. Second, Phase Estimation is a 
general framework to design and study quantum algorithms. Many known quantum algorithms 
that have exponential speedup over their classical counterparts can be described in terms of the 
phase estimation algorithm [201 Now our result implies that this is not an accident; actually 
all efficient quantum algorithms can be done in this framework. 

Measuring the energies of a quantum observable — or in a mathematical language, measuring 
eigenvalues of a Hamiltonian — is one of the most important tasks in quantum physics. The 
problems related to estimating/approximating eigenvalues of Hamiltonians or unitaries has been 
considered in the literature, including how to make the measurement for special cases [Tl [T?)l ITK] . 
or using the measurement to solve some computational tasks 30 , 8 . In this paper, we provide a 
precise complexity-theoretic framework for studying these problems, and show that the sampling 
versions of these two basic problems are actually BQP-complete. 

The remainder of this paper is organized as follows. Section |21 is devoted to basic notations 
used in the paper, and the definitions of the three problems we are studying in this paper. 
In Section |31 we prove that the three problems are BQP-hard. In Section ^ we show that 
the problems are in BQP. The paper concludes with Section |S1 which mentions some open 
problems for future work. 

2 Preliminaries and definitions 

A quantum register consisting of n qubits is described mathematically by a tensor product 
Hilbert space Ti = (C^)®". A state of the quantum register is given by a unit vector 
Transformations of the states are described by unitary operations acting on Ti.. Every uni- 
tary operation has to be composed out of elementary gates for implementation purposes. A 
measurement with respect to the orthogonal subspaces S'l, 5*2 (where 5*1 © 5'2 = H) causes the 
system to collapse to Pi | V-) / 1 1 Pi | V') 1 1 or P2 1 V-) / 1 1 A | V") 1 1 , with probabihty | ] Pi | V') | P and | ] P2 1 V') P , 
respectively, where Pi is the projector onto the subspace Si. For a comprehensive introduction 
to quantum computing, please refer to the textbooks (23 1^ . 

2.1 Definition of the problems 

The standard definition of Local Hamiltonian and Phase Estimation, two problems extensively 
studied in quantum complexity theory and in quantum algorithm engineering, are as follows. 

^Strictly speaking, approximation of the Jones polynomial is neither a language nor a decision problem. Here the 
concept of completeness is generalized in both their and our work: A computational task is BQP-complete if it can 
be solved by a polynomial time quantum computer, and any BQP problem can be reduced to this problem. 

■^Recently another problem (Consistency of Local Density Matrices) is showed to be also QMA-complete |24|. 
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Local Hamiltonian Minimal Eigenvalue (LH) 

Input: 1) Two numbers a and b such that a < h and the gap b — a = i}{l/poly{n)). 2) A 
Hamiltonian H = Hj operating on n qubits, where j ranges over a set of size polynomial in 
n, and each Hj operates on a constant number of qubits. It is promised that either X{H) < a 
or X{H) > b, where X{H) is the minimal eigenvalue of H. 

Output: if X{H) < a, and 1 if X{H) > b. 
Phase Estimation (PE) 

Input: A unitary matrix U, given by black-boxes of controUed-U', controlled-i/^^ , . . ., controUed- 
operations, and an eigenvector \u) of U with eigenvalue e^'^"^ with the value of (p G [0, 1) 
unknown. 

Output: An n-bit estimation of (p. 

The paper will study several average and sampling variants of them. For the sampling 
versions, we need a notion of the approximation of a probability distribution. Suppose there is 

a probability distribution p on the sample space X = {xi, .... xm}, on which a distance function 
d has been defined. Then a probability distribution on X is said to (e, 6)- approximate p if 

there is a decomposition qi = Ylf=i Hj s.t each qij > and Vj, '}2i:d(x,,x,)<e'i'i-j ^ (1 ~ ^)P3- 
Intuitively, this means that if we draw a sample according to q, then it will be e-close to Xj 
with probability at least (1 — 5)pj. Sometimes we also say that a sample drawn according to 
the probability distribution q is an estimation of Xj up to e with probability at least (1 — 5)pj. 

Comment. Note that because of the e-approximation (induced by the distance d), the standard 

distance measures between two probability distributions, such as total variance (or any p-norm 
generalization of it), Hellinger distance, or KuUback-Leibler distance, do not qualify for our 
purpose. For example, suppose X = {0.01,0.02,...,!}, p{x) = 0.1 if a; e {0.1, 0.2, 1}, and 
q{x) = 0.1 ii X G {0.09, 0.19, 0.99}. Then intuitively these two distributions are very close 
to each other. But the total variance between them is 5 X^^, \p{x) — q{x)\ = 1, the maximal 
possible distance. (In this example, the Hellinger distance is 1, and the Kullback-Leibler distance 
is infinity. All of these are not desired because they treat different samples in the sample space 
as totally different objects. But in many applications such as those in this paper, the sample 
space has a natural distance measure on it.) 

Now we are ready to define the sampling versions of LH and PE as follows. Note that they 
are sampling problems rather than the standard decision or searching problems: on every input, 
a sampling problem is required to output some value with some probability. 

Local Hamiltonian Eigenvalue Sampling (LHES) 

Input: 1) Hamiltonian H — Hj operating on n qubits, where j ranges over a set of size 
polynomial in n, and each Hj operates on a constant number of qubits. Suppose the eigenvalues 
and the corresponding eigenvectors of H are {(A^, \rjk))}k satisfying |Afc| < poly{n) for each k. 
2) An estimation precision e = U{l/poly{n)). 3) A sampling error probability S = ^l{l/poly{n)). 
4) A classical n-bit string b e {0, 1}". 

Output: An estimation of Afc up to e with probability at least (1 — ^)|(6|ryfe)|^. 
Phase Estimation Sampling (PES) 

Input: 1) ([/), the description of a 2" x 2" unitary matrix U (which is a poly{n)-size quan- 
tum circuit). Suppose the eigenvalues of U are {A^ = e'^'^'^'^^}j=i,...,2'<- (where (pj G [0,1) 
for each j), with the corresponding eigenvectors {|?7j)}j=i,...,2" • 2) An estimation precision 
e = Q{l/poly{nj). 3) A sampling error probability 6 = Q{l/poly{nj). 4) A classical n-bit string 

6e {0,1}". 

Output: An estimation of (pj up to e with probability at least (1 — ^)|(6|ryj)|^. 
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One could also consider to estimate the average phase instead of sampling the phase. One 
caveat here is that in the Phase Estimation and the PES problems, it is implicitly assumed that 
the estimation is in a circular sense, i.e. if the fu is very close to 1, then an answer very close to 
is considered correct.^ But when we consider the average phase, this becomes a problem: if we 
have two unitary matrices, one with all eigenvalues close to -1 and another with all eigenvalues 
close to 1, then we may not be able to distinguish them since the average phase of the second 
matrix could also be close to 1/2 due to the cancellation of those small and large phases. For 
this reason, it is more natural to consider the average eigenvalue of a unitary matrix. We define 
the Local Unitary matrix Average Eigenvalue (LUAE) as follows. 

Local Unitary matrix Average Eigenvalue (LUAE) 
Input: The same as that of PES. 

Output: An estimation of the average eigenvalue A = l(^l'7j)P'^j to precision e with 

probability at least 1 — 6. 

The main theorem of this paper can be written as 
Theorem 1 LHES, PES and LUAE are all BQF -complete. 

3 BQP-hardness of the problems 

In this section we will show that the three problems LHES, PES and LUAE are all BQP-hard. 

3.1 LHES is BQP-hard 

Theorem 2 BQP C BPP^^^^ 

Proof For any L e BQP, there is a uniform family of polynomial size quantum circuits with 
e-bounded error (for a small constant e) that decides if a: S L or a; ^ L as depicted in Figure Q 
For 71-bit inputs, denote by U the corresponding quantum circuit and suppose the size of U is 
M, which is bounded by a polynomial in n. Further suppose that the computation is described 
by U\x,0) = cex,o\0)\ipx,o) + ctx,i\l)\'4'x,i) , where is the initial state of the ancillary qubits, and 
\ipx.o) a-nd \ipx,i) are pure states. After the U transform, the first qubit is measured and the 
algorithm outputs the result based on the outcome of the measurement. As the requirement of 
correctness, we have |ax,oP < e if x G L, and |q!2;,iP < e if cc ^ L. 

2^ 



Figure 1: Circuit U 



^This can be seen from the analysis of the Phase Estimation algorithm; see, for example, )25l . 



— 
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We now construct a local Hamiltonian H encoding the circuit U and binary strings b encoding 
the inputs x such that eigenvalue sampling applied to H and b yields significantly different 
probability distributions for the two cases oi x € L and x ^ L. To this end, we construct the 
circuit V in Figure 13 where we add a qubit |r) to store a copy of the output of U. 



|0) 



u 



Figure 2: Circuit V 



Suppose the circuit is decomposed as U — Um • ■ ■ C^i, where each Uj an elementary gate. 

. , M, Vm+1 be 
1, the number 



Then [/t = • J/]^. Let Vj be the j-th gate in V, i.e. Vj = C/, for j 1, . 



,2M + 1. Let N = 2M 



the CNOT gate, and V, = C/2V+2-J j = M + 2, 
of gates in V. Attach a clock register \-)t to the system. Define the operator 



VN(E>\0)t{N-l\t. 



(1) 



Note that F is an 0(log A^)-local operator; we will remark how to slightly modify it to be a 
4-local operator at the end of the proof. Define 



\^,,,) ^ Fm0)\x,O)\0)t) 
for j > 0, where means that F is applied j times. Then for j = 0, 2N — 1, 

j=0,...,M 



(2) 



0)(g>Uj---Ui\x,O)(g,\j)t 

0) ® ■ ■ ■ uIiPqUm • • • t/i \x, 0) ® \j)t 

1) ® t/2V+2_^. • • • uIjPiUm ---Ui \x, 0) ® \j)t j = M + 1, . 

0) t^l-2M • • • uIPoUm • • • C/i \x, 0) ® |j - N)t 

1) ® t/]_2M • • • uIPiUm ---Ui \x, 0) (g) \j - N)t i = 2M + 1, 
0) ® t/2W-j • • • Ui\x, 0) ® \j ~ N)t i = 3M + 2, 



,2M 
. , 3M 



(3) 



where Pq and Pi are the projectors onto the subspaccs of the first output qubit of U being 
and 1, respectively. It is also easy to see that \ipx,j) — \Vx,j mod 2n) for j > 2N. 

Note that for different j and j' in {0, ...,27V — 1}, \(px,j) and \ fx,j') are orthogonal if \j — 
j'\ ^ N due to the clock register. Also note that Um---Ui\x,0) = U\x,0) = ax,o\0}\'4'x,o} + 
ax,i\^)\4'x,i) by definition. Therefore, we have {tp Wx,N+3) = {x,0\u1...uI^PoUm-Ui\x,Q) = 
|q^£!;,oP for i — 0, — 1. To summarize, we have 



{Vx,jWx.j') 



\Ol.xfi 



(4) 



Define the subspace Sx = span{\Lpx,j) '■ j = 0, ...,2A — 1}. The key point here is that 
though F has an exponentially large dimension, F\s^ is of only polynomial dimension. Now 
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if cnxfi = Ij then Sx has dimension of N, and F has a period of N on Sx- ActuaUy, F\s^ 



is just a shift operator on the basis {\'Px,j)}j=o,...N-i, 



F\(pxj 



\VxJ + l mod n) 



It is 



not hard to see that this operator has eigenvalues = with corresponding eigenvectors 



l^j.) = ^^X]j=o '^N ''\'Px,j)j where ujn = e^'^^l^ . Also note the fact that 

|(¥'.,o|a>l' = l/^^ 



(5) 



for any k — 0, TV — 1. 

Now we consider the general case of ax.o 7^ 1, in which case Sx has dimension 2iV and F has 
period of IN . Note that if x ^ L, then \ipx.N+j) is very close to l^fixj), if a; G L, then \ipx^N+j) 
is almost orthogonal to \<fx,j)- See Figure |21 for an illustration. 



l<Px,N> 

^l9x,0> 



l9x,N+l> 

^l9x,l> 
(a) X ^ L 



l9x,2N-l> 
*l9x,N-l> 



x,N ■ 



,l9x,N+l- 



,l9x,2N-l> 



l9x,0> 



l9x,l> 

(b) X e L 



l9x,N-l> 



Figure 3: The 2N vectors in the two cases 



To find the eigenvalues and eigenvectors of Fj^^, define 



Wlfx.j) + \'^x,N+j)\\ 



Wx,. 



Wx.N+j) 



(6) 



Wx,o) - \ipx.N+j)\\ 

for j = Q,...,N — 1. Then first, because {(pxj\^x,N+j) = |q^x,oP is a real number, we have 
{4>x,jhx,j) = 0. Together with Equality Q, we know that {\(f>x.,o), \4'x,n-i), hx.,o), hx.N^i)} 
form an orthonormal basis of S^.. Second, since 



l^fixfi) ^ \Vx,l) ^ 



F F 
— > \'fx,2N-l) — > WxS)) 



we observe that 
and 



F F 

\(l>xfi) — > \4>x.i) 



\lxfi) |7x,l> ^ 



F F 

\4'x,2N-l) — ^ l^'x.o) 



F . . F 

— > nx,2N-i) — > 



llxfi) 



Therefore F\s^ = Fi ® F2, where Fi and F2 act on Sx,+ — span{\(f>xfl) , ■■■,\(I)x,n-i)} and 
Sx,- — span{\jx.o)T ■■■Alx,N-i)}, respectively, with the matrix representations (in the basis 
{|</>a;j)} and {|7a;j)}, respectively) as follows: 



/O ... 1\ 
1 ... 
1 ... 



\0 



F2 



1 oJ 



/O 
1 
1 



Vo 



-1\ 





1 0/ 



(7) 
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So Fi is just a shift operator we mentioned before, which has eigenvalues /ife = and eigen- 
vectors \rik) = -^J2j'=o^'^N'^''\'f'x.j)- It is also not hard to find that F2 has the eigenvalues 

j/fc = L^'^^^'^ with the eigenvectors — -^J2f=o^ ''^N^'^^^^^^''\"fx.j)- P+ and P_ be the 
projectors onto 5*+ and S^, respectively. Note that 

IIP I M|2 _ II \ - \y^x.N) ||2 _ 2 - {yyx,o\y^x,N) - (y^x.jvlt/^a.o) _ 1 " |«:r,Op 



Thus ||P+|(^j;,o)|P = 1 - ||P-|(/9:e,o)||^ = (1 + |aa;,oP)/2. By this, we can derive that 

|(¥'.,o|?7fe)l'-^(l + |a..on, l(V'.,okfc)|' = ^|a.,i|'. (9) 

for any k = 0, N — 1. 

Now the local Hamiltonian is constructed a,s H — F + . It is easy to verify that _ff is a 
local Hamiltonian. And further, suppose the eigenvalues of F are {k^} with the eigenvectors 
ll'i/'j)}, then the eigenvalues of H are just {kj + n*} with the same corresponding eigenvectors. 
Now the BPP algorithm with the LHES oracle is as follows. 

On input x, 

1. Feed (H, l/4:N, 1/100, OxOO) as input to the LHES oracle, getting an output a 

2. If \a\ > 1, 

go back to StepQ 

3. If minfc=o,...,JV-i \a - 2 cos ^'^('^+^2) | ^ niinfc=o,...,Af-i |a- 2cos2^| 

output "a; G L" 
else 

output "a; ^ L" 

The algorithm is correct with probability at least i(l — (5)(1 — e) for x E L, and correct with 
probability at least (1 — 5)(1 — |) for x ^ L . Actually, in Step^ the oracle will output an 
estimation of 2cos^^ (up to an additive j^) with probability at least 27v(-'^ ~ + l«2;,oP), 
and output an estimation of 2 cos 2^(^+^/2) (^p to an additive ^) with probability at least 
2^(1 — (5)|Q;a;,iP, for any fc = 0, — 1. Note that now the case of ax.o = ^ can be consistently 
viewed as a special case of ax,o 7^ 1 by comparing Equality jSJl and (Q. 

li k G [A/6, A^/3]U[2A^/3, 5N/6], which happens with probability 1/3 for a uniformly random 
k € {0, A^ — 1}, then the algorithm will proceed to the outer "else" branch in Step|21 Note 
that for 0,6'' e [tt/S, 27r/3]U[47r/3, 57r/3] and \0-0'\ = 27^, we have 1 2 cos 6* -2 cos 6*' | > ^ > 2^- 
So if X G L, then with at least probability ^(1 — (5)|aa;_ip > ^(1 — S){1 — e), we can get a value 
" < 47v"'^^°^'^" 2 '^'^^ 2TT{k+i/2) some fe, but "> ^-far" from any 2 cos so the last step 
will catch this and output the correct answer. 

The case oi x ^ L can be similarly analyzed. Also, it is easy to see that the expected 
running time of the algorithm is polynomial in the input size, which completes the proof for the 
0(log n)-local LHES reduction. 

Finally, to obtain a 4-local LHES, we replace \i) by \ei) = |0...010...0) for « = . . . , A^ - 1, 
where the only 1 appears at coordinate i. Modify the operator F to be 

JV-l 

F=^Vj® |ej)(ej_i| + Vn® |eo)(eAr_i| . (10) 
Note that |ej)(ej_i| and |eo)(eAr_i| are 2-local. The remaining proof passes through. □ 
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3.2 PES and LUAE are BQP-hard 
Theorem 3 BQP C BPP^^^, BQP C BPP^ 



Proof For any L G BQP, there is a uniform family of polynomial size quantum circuits with 
e-bounded error (for a small constant e) that decides if a; G L or a; ^ L as in Figure ^ For any 
input X G {0, 1}", suppose the corresponding circuit U operates on iV = poly{n) qubits, taking 
X and some ancillary bits — 0^^" as the input. After the circuit computation, measuring 
the first qubit of the output [/|a;)|0) give the correct result with probability great than 1 — e. 
Suppose U\x)\Q) = Q!j;,o|0)|?/'2;,o) + aa:,i\l)\il)x,i) , then |a:E,oP < e if a; G i, and \axs? < e if 
X ^ L. 

Now we construct another circuit V as in Figure 01 where Z is the Pauli-Z matrix 



^2) 



|0) 



U 



z 



[/t 



Figure 4: Circuit V 



It is easy to see that 

y|x,0) =C/^K,o|0)|7/^,^o) -«.,i|l>IV'.,i)) (11) 
|x,0)-2a,,iC/t|l)|^,^i) 
\-|x,0) + 2a,,oC/^|0)|^,,o) 

Therefore, |x, 0) is almost an eigenvector of the eigenvalue —1 or 1, depending on whether x (z L 
or not. So we can use the PES or LUAE oracle to distinguish between these two cases. 

1. BQP C BPP^^-^: We run the PES oracle on input {{V), 1/8, 1/100, xO^""), and the oracle 
outputs a (sampled) value 9. Accept if 1/4 < < 3/4, and reject otherwise. 

To see why this is correct, write V as = J2j ^j\Vj){'']j\^ where Xj = e^'^*'^^ are eigenvalues 
and \r]j) are the corresponding eigenvectors. Consider x <^ L first, and the other case x € L 
is symmetric. Denote xO^^" by b and decompose it as \b) = Then V\b) — \b) = 

- ^)Pj\Vj)- By this implies 

^(A, - l)/3,h,) = -2a,.iC/t|l)|^,, 1). (13) 
j 
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Since the PES oracle outputs an estimation of Lpj up to 1/8 with probability at least (1 — 
i5)|(6|77j)p with 5 = 1/100, the success probability of our algorithm is 



Pr[6'< l/4or6'>3/4]> ^ {1 - S)\{b\r]j)f 

j:ipj<l/8 or i^j>7/8 

i:|Aj-l|<2sin(7r/8) 



= {1-5 
= {l-S 

>{l-5 

>il-6 
= {1-S 
= {1~S 

= {1-6 
= {1-6 



j:|Aj-l|>2sin(7r/8) 



E 



j:|Aj-l|>2sin(ir/8) 



4sin2(V8)"^^' ' 



(i-ET^S7^i/3.n 



4sin^(7r/8) 



(1 



4sin^(7r/8) 
1 

4sin2(7r/8) 

sin2(7r/8) 
e 



|2a,,it/^|l)|V,,i)f) 



sin^(7r/8) 



(14) 
(15) 
(16) 

(17) 

(18) 
(19) 
(20) 

(21) 
(22) 



where (|20|) is because of equality (|13|l . Therefore the error probability is less than 0{e). The 
same arguments can show the correctness for the case of a; e L, which completes our proof for 
the first part. 



2. BQP C BPP 



LUAE. 



For the circuit V in Figure 0] it is easy to show that A is close to 



^1 (and 1) if a; G L (and x ^ L). First observe that for any b, we have A = {b\U\b), which 
can be shown by writing both U and b in terms of |»yj)'s. Now for b = xO, ii x ^ L, then 
|A- 1| = \{b\U\b) -1\ = \{b\U\b) - {b\b)\ <p\b) - \b)\ = 2\a^,i\ < 2^. Similarly if a; € L, then 
|A + 1| < Thus a good estimation to A suffices to distinguish between x & L and x ^ L. D 



4 The problems are in BQP 

In this section, we will prove that LHES, PES and LUAE are all in BQP. We will ffi'st review 
the standard algorithm for Phase Estimation, then observe that the same algorithm actually 
gives the desired PES solution. We then use it to show an algorithm for LHES. Finally, we 
prove that LUAE is in BQP. 

4.1 Phase Estimation and an efficient quantum algorithm for PES 

Phase Estimation can be solved by a quantum algorithm as follows (see, for example, |2S1)- The 
working space has two registers. The first register consists ol t = n + [log(2 + 1/2^)] qubits 
and is prepared in |0. . .0). The second register contains the eigenvector \u). Measuring (p in 
the first register after carrying out the transformations described below gives the desired n-bit 
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estimation of ip with probability of at least 1 — 5. 

|0)«*|«) (23) 
1 

// apply the Fourier transform (24) 

V 2* ■ n 

/= 5Z l-^)^"*!^) // ^PPly t-^^*^ controlled powers of U (25) 



2* .^0 



2*-l 

s'-^'^Ij)!^^) (26) 



-^\(p)\u) 1 1 apply the inverse Fourier transform (27) 

The following observation says that the same algorithm actually works for PES. 

Fact. If we feed |0) |6) instead of |0) |u) as input to the above algorithm for the Phase Estimation 
problem and let t = [log + [log(2 + , then the measurement of the first register gives the 
desired sampling output. This implies that PES can be solved by a BQP machine. 

This actually holds not only for \h) but also for a general state \r]). To see why this is true, 
write \ri) as X]j=i '^jl'7i)' then by the linearity of the operations, the final state is aj\ipj)\rij) . For 
more details, we refer the readers to |2S] (Chapter 5). Note that to implement the controUed- 
U"^' operations for j = 0, 2* — 1 in the above algorithm, we need to run U for 2* times, which 
can be done efficiently since t = [log + [log(2 + ^)] and e = / poly{n)) , 5 — / poly{n)) . 

4.2 LHES in BQP 

Theorem 4 LHES can he implemented by a uniform family of quantum circuits of polynomial 
size. 

Proof By a simple scaling {H' = H/ A = J^j Hj/ A where A — max^ |Afc| = poly{n)), we can 
assume that all the eigenvalues of H satisfy |Afc| < 1/4. The basic idea to design the BQP 
algorithm is to use Phase Estimation Sampling on e'^'^^^ . Note that e^^'^ is unitary, and if 
the eigenvalues and eigenvectors of H are {Afc, \r]k)}, then those of e^'^'^ are just {e^''*'^*, \rik)}. 
Therefore, it seems that it is enough to run the PES algorithm on (e^^*^, e, 5, 6), and if we 
get some A > 1/2, then output A — 1. However, note that H is of exponential dimension, so 
g27riff j-^Q^ ready to compute in the straightforward way. Fortunately, this issue is well studied 
in the quantum simulation algorithms, and the standard approach is the following asymptotic 
approximation by the Trotter formula [281 1101 125j or its variants. Here using the simulation 
technique, we obtain 

/ 2^,:i:^.ff,/m\'" ^ j-Qg2^,ff,/m 1 +0{l/m) (28) 



Now we run PES on {b, e, S/2, e^'^^^). Whenever we need to call e^'^*^, we use Y\j e^'^*^^/™ for 
m times instead. Note that such substitution yields 0(l/m) deviation, sot — log ^ + 0(1) calls 
yield 0(—^) < — ^ deviation for some constant c. Let rn = thus the final error probability 
is less than | + = S, achieving the desired estimation and sampling precisions. □ 

Comment. From the proof we can see that as long as e^'^*^ can be simulated efficiently from 
the description of _ff , we can sample the eigenvalues of H as desired. Since sparse Hamiltonians, 
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which contains local Hamiltonians as special cases, can be simulated efficiently we know 

that if we modify the definition of LHES by allowing H to be sparse, then it is also BQP- 
complete. 

4.3 Local unitary matrices average eigenvalue estimation 
Theorem 5 LUAE is m BQP. 

We can use the BQP algorithm for PES to get 0{l/{e5)) (independent) samples '■P{i), V(m) 
and use the sample mean ^ ~ Sj=i e^'^*'''*^' as an desired estimation of A, which can be 
proved by studying the Phase Estimation algorithm in more details But here we shall give 
a different algorithm whose analysis is much simpler. Recall that A = {h\U\b). So after we 
apply the circuit U to the problem becomes to estimate the inner product of two quantum 
states \b) and U\h), which can be done by the standard SWAP test. For example, Yao gave 
an estimate of up to e with error probability i5 by applying the SWAP test in [7] on 

0(^log|) pairs of \v)). Here for the special two states \b) and U\h), it turns out that we 
can slightly decrease the cost to log j-), by using the following lemma proved in |31j: 

Lemma 6 Let U he a quantum circuit of length 0{poly{n)) acting on n qubits, and let lip) be 
a pure state of the n qubits which can he prepared in time 0{poly{n)). It is then possible to 
sample from random variables — for which 

E[X + ir] = (^|i7|i/.) (29) 

in 0{poly{n)) time. 

Proof (of Theorem O Run the sampling algorithm in the lemma to get m = 0(^log|) 
samples (Xi, Yi), (X„, y„). Let X = ^ LJli and Y = ^LJli^j- Use the sample 
mean A = X + as an desired estimation of A. Noting Var[X] < 4 and Var[y] < 4 (by the 
definition of variance), we have by Chernoff's bound that 

Pr[|A - A| > e] < Vy[\X - Re{X)\ > e/2] + Pr[\Y - Im(X)\ > e/2] < S. (30) 

as desired. □ 



5 Discussions and open problems 

One might wonder why not to consider the average eigenvalue of a local Hamiltonian. Actually, 
since {b\H\b) — {b\ J2j Hj\b) — and each {b\Hj\b) can be easily computed even deter- 

ministically, we can obtain the exact average eigenvalue of a local Hamiltonian detcrministically 
in polynomial time. 

However, as shown very recently by Janzing and Wocjan, if we generalize the problem to 
estimating the average eigenvalue of -ff™ where m = polyin) is part of the input, then the 
problem of estimating {b\H"^\b) is BQP-complete |17) . 

We can also consider the complexity of unguided version of the problems studied so far. That 
is, 6 is not part of the input, and we want to sample the eigenvalue/phase with equal probability 
or estimate the average eigenvalue under the uniform distribution over all the 2" eigenvalues. We 
use LHESu, PESu, LUAEu to denote the corresponding problems. We basically do not know 
anything about these problems yet, except some trivial facts like LUAEu G BQP due to the 
following simple observation: Eb[(6|C/|6)] = ^ Eb{b\U\b) = ^ J^bj \{bh)\^^j = w = A. 

(We can just uniformly choose b e {0, 1}" and run the algorithm in Lemma|Sl) 

''Directly applying the definition of the PES problem as a black-box is not enough: the mean of the output of the 
PES algorithm may be far away from A because of the exponentially many eigenvalues and for each eigenvalue the 
output probability is e away from the correct one. 
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